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We provide a rigorous generalization of the quantum adiabatic theorem for open systems described 
by a Markovian master equation with time-dependent Liouvillian C{t). We focus on the finite system 
case relevant for adiabatic quantum computing and quantum annealing. Adiabaticity is defined in 
terms of closeness to the instantaneous steady state. While the general result is conceptually similar 
to the closed system case, there are important differences. Namely, a system initialized in the 
zero-eigenvalue elgenspace of C{t) will remain in this eigenspace with a deviation that is inversely 
proportional to the total evolution time T. In the case of a finite number of level crossings the scaling 
becomes with an exponent rj that we relate to the rate of the gap closing. For master equations 
that describe relaxation to thermal equilibrium, we show that the evolution time T should be long 
compared to the corresponding minimum inverse gap squared of £-{t). Our results are illustrated 
with several examples. 


I. INTRODUCTION 


The origins of the celebrated quantum adiabatic the¬ 
orem (QAT) date back to Einstein’s “Adiabatenhy- 
pothese”: “If a system be affected in a reversible adia¬ 
batic way, allowed motions are transformed into allowed 
motions” [1]. Ehrenfest was the first to appreciate the 
importance of adiabatic invariance, guessing—before the 
advent of a complete quantum theory— that quantum 
laws would only allow motions which are invariant under 
adiabatic perturbations [2]. The more familiar, modern 
version of the QAT was put forth by Born and Fock al¬ 
ready in 1928 for the case of discrete spectra [3]. Since 
then a series of increasingly sophisticated techniques have 
been developed in order to generalize the QAT to include 
degeneracy, unbounded models, continuous spectra, and 
exponential error estimates [4-11]. 

This long history of adiabatic theorems is almost ex¬ 
clusively concerned with closed systems undergoing uni¬ 
tary evolution. Previous approaches to formulating an 
adiabatic condition for open quantum systems [12] have 
focused on a Jordan block decomposition of the dissi¬ 
pative generator [13], the weak coupling limit [14, 15], 
zero temperature [16], or on a noiseless subsystem de¬ 
composition [17, 18]. Here we prove that in analogy to 
the closed system case, where the system follows the in¬ 
stantaneous (pure) eigenstates of the Hamiltonian in the 
adiabatic limit of arbitrarily large total evolution time T, 
the open system follows the instantaneous steady state 
(ISS) of the Liouvillian. In doing so, we extend the sem¬ 
inal closed system result by Kato [4[ to infinite order in 
IjT. 

Rigorous extensions of the adiabatic theorem for gener¬ 
ators of contractive semigroups, similar to ours, have also 
appeared in the mathematical literature [19-22]. Our fo¬ 
cus is on estimating the adiabatic error in terms of the 
physical parameters of the theory, thus making the re¬ 


sult more suitable for applications. We demonstrate that 
in the case of thermal baths satisfying the Kubo-Martin- 
Schwinger (KMS) condition [23], where the ISS is the 
instantaneous Gibbs state, a sufficient criterion for adi¬ 
abaticity is T ^ where Amin is the smallest Li- 

ouvillian gap in absolute value. Our QAT also allows 
for a finite number of level crossings in the Liouvillian 
spectrum, for which we demonstrate that the error scales 
asymptotically as T~^ with a known exponent r] € (0,1) 
related to the rate at which the gap closes. This setting 
is directly relevant to recent theoretical and experimental 
work on quantum annealing [15, 24-26], where the Liou¬ 
villian gap may close at the end of the evolution, and we 
verify this prediction using numerical simulations. 


II. INSTANTANEOUS STEADY STATES 

We assume that the evolution of a d-dimensional (d < 
oo) quantum system with state p(t) can be described 
by a linear, time-local master equation dp/dt = CT{t)p, 
where T is the total evolution time. We also assume that 
Ct{sT) = £(s) where s = t/T G [0,1] is a rescaled, di¬ 
mensionless time coordinate, and C{s) is T-independent. 
Setting pt{s) = £t{s,0)pt{0), the evolution operator 
^r(s,so) satisfies 


£^{s, So) = T£{s)£t{s, So) 


( 1 ) 


with f(s,s) = 1 (we drop the subscript T from now 
on; we also write £{s) for £i(s,0) for simplicity), and 
where the prime denotes dg ■ We are interested in 
the solutions of Eq. (1) for large T. We further 
assume that the Liouvillian C{s) can be written in 
Lindblad form for all s, i.e., C{s)9 = —i[dd(s),*] -I- 

Liis) • lI{s) - ^{lI{s)Li{s),»} , where id( s) is the 
system Hamiltonian and {Li{s)} are the Lindblad oper- 
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ators. Equation (1) then describes a Markovian master 
equation with a time-dependent Lindblad generator, and 
the corresponding evolntion operator £(s 2 ,si) is a com¬ 
pletely positive trace preserving (CPTP) map for any 
S 2 > Si [15, 27, 28]. We formulate our results in terms 
of Lindblad operators and CPTP maps, bnt in fact all 
our results are valid in the more general case where C{s) 
generates a contractive semigroup (i.e., |j ||< 1 Vs 

and t > 0). A special role is played by the ISSs, i.e., the 
states in the kernel of C{s). In the time-independent case 
[£(s) = £, Vs], it follows from the CPTP property that 
any initial state evolves to Ker£ in the long time limit 
]29]. Let ns denote by P{s) the (instantaneons) spectral 
projection of £(s) with eigenvalne zero. The Lindblad 
form guarantees that zero is a semi-simple, possibly de¬ 
generate, eigenvalue of £(s) (see Appendix A for a proof), 
and so there are no nilpotent terms in the zero sector, i.e., 
Cis)P{s) = P{s)C{s) =0. 

We are now ready to informally state the QAT for open 
systems: If a system is initialized at s = 0 in Ker£(0), 
the final state at s = 1 will be close to Ker£(l), provided 
the Lindbladian changes sufficiently slowly. 

In principle one conld formulate an open-system QAT 
considering other (non-zero) eigenvalues of £(s). How¬ 
ever their corresponding invariant snbspaces contain no 
physical states, so that the physical interest in snch a 
generalization is questionable [30]. We proceed to rig- 
oronsly establish the QAT and identify the timescales it 
entails. 


III. GAPPED CASE 

We start by assuming that the zero eigenvalue is sep¬ 
arated by a finite gap Amin from the rest of the spec¬ 
trum a{C{s)) for all s S [0, 1], i.e., dist ](t(£(s))\{0}, 0] = 
A(s) > Amin > 0, a condition we relax later. The ideal 
adiabatic evolution is represented by an operator V (s) 
that satisfies the intertwining property: V{s)P{0) = 
P{s)V{s). It is well known [4] (see Appendix. B) that 
a possible choice for P(s) is given by the solution of 
the differential equation W(s) = ]P'(s), P(s)]P(s) with 
P(0) = 1. We are interested in quantifying the devi¬ 
ation of the actual evolution, governed by the CPTP 
map £(s), from the ideal adiabatic evolution. However, 
V(s) is not CPTP in general. Instead, we can prove that 
W{s) := P(s)P(0) is a CPTP map, since it can be writ¬ 
ten as a product of projectors (CPTP maps): W{s) = 
limAr_>oo P(s) • • ■ P{2s/N)P{s/N)P{0) (see Appendix C) 
or Proposition 3 in [22]). Therefore, to state the QAT we 
wish to bound the deviation from the ideal adiabatic evo¬ 
lntion projected to Ker£, ||£(s)P(0) — P(s)P(0)||, in the 
large T limit [31]. 

To proceed we introduce the reduced resolvent S{s) = 
lim^_>o Q(s)(£(s) - z)~'^Q{s), where Q{s) = 1 - P(s) 
[32]. We fnrther assnme that C is m times differentiable, 
and let A„_|_i(s) = S'(s)X(,(s), with Ai(s) = S{s). Un¬ 
der the additional simplifying assnmption that the ISS is 


unique we then prove the following in Appendix D using 
integration by parts: 

m 

]£(s)-U(s)]P(0) = ^^ (2) 

n—1 

~ ^ lo ’ 

£l„ = £{s,a)Xri{a)W'{a)\l- f da£{s,a)Xn{a)W''{a) . 

Jo 

The general resnlt, valid also for degenerate kernels, is 
given in Appendix F. 

It tnrns ont that Eq. (2) is valid in the to = 1 case 
even withont requiring that the ISS be unique, and as we 
show in Appendix D we can bound the deviation from 
the ideal adiabatic evolntion in general as 

\\[£{s)-V{s)]Pm<C/T . (3) 

A similar result has been derived in [22] where, however, 
the constant C is left undetermined. We show that the 
constant C, independent of P, can be taken to be 

C=||P(s)||||P'(s)|i + ||5(0)||||P'(0)|| 

+ sup \\[S'P' + SP''Ka)\\ . (4) 

Below we discnss how C relates to the physical parame¬ 
ters of the model associated with £. Inequality (o) states 
the QAT for open systems and implies the QAT in the 
standard form for states. 

To see the latter let ns initialize the system in a state 
p(0) in Ker£(0), i.e., p(0) = P(0)p(0). Then p{s) := 
U(s)/5(0) is an instantaneons steady state at time s, since 
Cis)pis) = Cis)Vis)P{0)m = Cis)P{s)V{s)piO) = 
0, i.e., p{s) € Ker£(s). Under the actual evolution 
the state is mapped to p{s) = £{s)p{0) and one has 

lip(s)-p(s)||i < |1 ]£(s)-U(s)]P(0)|||1p(0)|U < C/T. 

Namely, if T > C/e (e > 0) then the system is gnaran- 
teed to find itself e-close in norm to the instantaneons 
steady state, at the end of the evolution. 

IV. THE CLOSED SYSTEM LIMIT 

It is nsefnl to comment on how our result relates to 
adiabatic theorems for closed systems, described by a 
system Hamiltonian i7(s) with eigenvalnes En. 

First, if one is interested in initial states belonging to 
the —iEo level, one may simply set £(s) = —i[II{s) — Eo], 
as our formalism encompasses (with minor modifications) 
the case where £ is anti-Hermitian. In this way one re¬ 
covers the standard adiabatic theorem for closed system. 
The relevant gap is given by the eigenvalue closest to Eq 
in modnlns, i.e., \Ei — Eq\. The bonnd we obtained for 
the constant C in this case is similar to that given in 
Ref. [10], at least for what concerns the dependence on 
the gap. 
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Another possibility is to write £(s) = JC{s) = 

—The eigenvalues of /C are {—i{En — Em)}- 
There is a > d-fold degenerate zero eigenvalue arising 
from En = Em constituting KerA. The relevant en¬ 
ergy scale is determined by the next eigenvalue which 
is closest to zero in modulus. This is given by the small¬ 
est difference |if„ — Em\ with 7 ^ Em (non-zero since 
we assume d < 00 and hence discrete spectra). This is 
consistent with the previous result because in Ker/C one 
has the freedom to pick any state \n){n\ leading to a gap 
min^ \En — Em\- In this manner one obtains an adiabatic 
theorem for closed systems in the Liouvillian (superop¬ 
erator) formalism. We discuss the closed system limit 
further in Appendix E. 


V. THERMAL BATH 

It turns out that the open system version of the QAT 
can have additional structure that is absent in the closed 
system case. To demonstrate this we consider the im¬ 
portant class of Lindbladians generated by the interac¬ 
tion of a system with a thermal bath, for which we can 
make the bounds above more specific. As a result of the 
KMS condition such Lindbladians satisfy the quantum 
detailed balance condition [33, 34]. This fact has im¬ 
portant consequences, namely, (i) the Gibbs state is an 
ISS, i.e., £(s)pg(s) = 0 with pg{s) = exp[—PH{s)]/Z, 
where Z = Trexp [—(3H{s)] is the partition function; (ii) 
the generator C{s) is normal. 

Let us now show how we can relate C to standard quan¬ 
tities such as the gap and H'{s) using the assumption of 
a thermal bath. Assume for simplicity that pg(s) is the 
unique state in Ker£(s). For such thermal baths the pro¬ 
jector onto the ISS manifold is P{s) = |pg(s))( 1| Then 

p(")(s) = |p[?)(s))(l|, so that ||P("Hs)ll = IIPg^(s)||i- 
Thus, if H(s) is bounded with bounded derivatives, 
||P*'"^(s)|| is bounded for all s, and hence P' and P" 
in the constant C [Eq. ( )] do not introduce any sin¬ 
gularities. In addition, since C is normal, ||S'|| = c/A 
where the constant c depends only on the norm used 
[35]. Moreover, the identity S' = S'^CP + PC'— SCS 
(Appendix G) implies that ||S''|j < 3||S'||^||£'||. Thus, 
from Eq. (i), for thermal baths the dependence on the 
Liouvillian gap is C = 0(A~?^). Note that in the ab¬ 
sence of the KMS condition one only has C = 0(A“?,), 
implying T = 0(A“?^) as a criterion for adiabaticity 
in accordance with the closed system result of Ref. [10] . 
However, in particular cases the dependence can be even 
milder. For example, if C{s) is a unitary family, i.e., 
C{s) = e®^£(0)e“®^, with /C an anti-hermitian superop¬ 
erator, one has (Appendix H) P' = KP — PIC, so neither 
P' (nor P") depends on Amin- Moreover S' = ICS — SK. 
so that in this case C = 0(A“Ij,), as also shown and 
extensively exploited in [36, 37]. 

In addition, for thermal baths the constants appearing 
in Eq. ('*) bear an explicit dependence on H(s). E.g., 


one can show that ||pG^(s)||i < 2 /?A/([iJ'(s)] 2 )( 3 , where 
(•)g = Tr[pG*] is the thermal average (see Appendix I). 
This fact has important consequences for adiabatic quan¬ 
tum computation where the complexity of a computation 
is encoded into H{s) and depends on the system size L. 
In general we expect |1-P^"'H'S)II to display a stronger di¬ 
vergence with L, for some n, e.g., at (positive tempera¬ 
ture) phase transition points of H{s) [38] . 

When the gap Amin is very small and is attained in¬ 
side the interval [0, s], the constant C is dominated by 
the third term in Eq. (1), i.e., C ~ sup^ ||S''(CT)P'(cr)||. 
Using the above estimates for S' and P' we obtain 
c < 6c2/3||/:'(cr)Umax a/([-H'(ct)]2)g, maxA“f„, where the 
subscript “max” means that the corresponding quantities 
must be maximized over a € [0, s]. In other words, taking 
T > c^/3||/:'(cr)||max\/([-ff'(o-)]^)G.maxA“f,,/e guarantees 
adiabaticity up to an error 0(e) in trace norm. 

VI. CASE OF LEVEL CROSSINGS 

The gapped case is typical since a random Lindbla- 
dian will have a gap above zero with probability one for 
all values of s. However, symmetries may give rise to de¬ 
generacies, and so we would like to extend our result and 
consider the case where a finite number of level crossings 
with the zero eigenvalue may take place along the path 

[39[; 

Since singularities are only algebraic in the finite di¬ 
mensional case, it is reasonable to expect that, in the case 
of level crossing, one has || [£(s) — U(s)]P( 0 )|| ^ l/T’^ 
for large T, with a positive exponent rj < 1. We are 
interested in estimating rj for large T. For definiteness 
assume that at the level crossings the gap vanishes as 
Amin(s) — Vi{s — s*)“* with some positive exponents 
The analysis is detailed in Appendix J. The final result 
is 

^ D 1 

II [£(s) - U(s)] P(0)|| < E ^ ’ (5) 

“ J 1 + ai 

1—1 

where Di are positive constants. Clearly the asymptotic 
behavior of the right hand side is dictated by the smallest 
exponent rji, i.e., by the largest ai, and hence the most 
divergent of the N gaps. 

VH. EXAMPLES 

We now illustrate our results with a few examples. 

A. Example 1 

Let us first consider a time-dependent generalization of 
the amplitude damping master equation. The Lindbla- 
dian is £(s) = lC{s)+Co with lC{s) = —i[H{s), •], H{s) = 
m{s) ■ cr and Cq = 2'y [a~ • cr+ — (1/2) {(T+cr“, •}] (tT“ 
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Figure 1. QAT for the gapped case using the model of Exam¬ 
ple 1. (a) Absolute values of the eigenvalues of the Lindbla- 
dian (zero is not shown). The eigenvalues show square-root 
singularities at s ~ 0.88 and s ~ 0.94. (b) Trace norm dis¬ 
tance between the actual state and the ISS for increasing T. 
The blue line is given by ||p(l) — p(l)|| ~ 7.88/T. Parame¬ 
ters: mx{s) = 1 — s, my{s) = 0, mz{s) = s/150, 7 = 1/2 arb. 
units Initial condition is p(0) as given in Eq. ( 6 ). 




Figure 2. Illustration of the QAT for a single qubit coupled to 
a thermal bath (Example 2). (a) Gapped case with A = gcr^ ■ 
The at (blue line) gives ||p(l) - p(l)|ji = 148.5/r° ®®®°. (b) 
Gapless case with A = . Now the ht gives ||p(l) —p(l)||i = 

1.910/r'’'^^'*. Fits are for T > 10®'"^. Parameters: g = 10“^, 
Ljj; = LOz = ~l/2 arb. units, /3 = 1 arb. units, 7(0;) = 


Pauli matrices and cr^ = cr^ ± ia'^). The steady state 
manifold is one-dimensional. The ISS is given by the 
solution of C{s)p{s) = 0 and is (in the tr^ basis) 

- _ 1 / — ml —m-{2mz + iy) \ 

^ c \ —TO+( 2 m 2 — 17 ) + 3ml + 7 ^ y ’ 

where c = 2(m^ -|- ml) + 7 ^ and m± = mx ± imy. A plot 
of the absolute values of the Lindbladian’s eigenvalues is 
shown in Fig. 1(a). Parameters are chosen to illustrate 
a phenomenon which is not possible in the unitary case. 
Namely, the eigenvalues can have algebraic singnlarities 
of the form (s—s*)’" with r a non-integer rational number. 
Also shown, in Fig. 1(b), is the expected scaling of ||p(l) — 
p(l)||i asC/T. 


B. Example 2 


Consider a qubit with system Hamiltonian H{s) = 
m{s) ■ (T interacting with a heat bath at inverse tem- 
peratnre /3, so that the total Hamiltonian is Htot(i) = 
H{t) + Hint + Hq, with Hint = A® B, where A (B) is a 
system (bath) operator, and Hb is the bath Hamiltonian. 
For H(s) of the form of Landau-Zener driving, exact ex¬ 
pressions for the transition probabilities for a particu¬ 
lar zero-temperature bosonic bath have been obtained in 
Ref. [40]. To treat the more general case we nse the time- 
dependent Lindblad master equation approximation for 
a slowly varying system Hamiltonian [15]: £ tot(^) — 

/C(t) + £(t), where /C(t) = -i [H{t) + i?Ls(0>*]- 

The dissipative part, specialized to the single qnbit 
case, reads 


£{t) = J27{uj)\AUt)>AUt)^ - I {A^(t)tA^(t),.} 


(7) 

with rates 7 ( 0 ;) = The 

Lindblad operators are given by the Fourier reso¬ 
lution e^rH(t) ^ Ay,{t), where w 


are the Bohr frequencies of H{t). H^sit) = 
S{uj)Ay,{t)^A^{t) is the Lamb shift Hamiltonian, 

with S{uj) = duj''f{uj')V ^ ^, where V denotes 

the principal valne. Moreover, 7 q(— cn) = 
as a conseqnence of the KMS condition. This implies 
that the generator C{t) together with IC{t) and /Cls]^) = 
—i [H-Bsit), •] all commnte at the same time. 

The corresponding Lindbladian has the following in- 
stantaneons eigenvalues (derived in Appendix K): 

A = { 0 , -7(5) lAoij' (1 + e-^^), -F ± ip] , (8) 

with 2r = 7(0)(|Aoo|^ + [Anj^) -b 7(5) |Aoi|^ (1 -b e~^^) 

and p = 5- 5(0)(|Aoo|' - [An|') + [Aqi[ 5(5) - 5(-5)], 
where 5 = 2||m|| is the instantaneons Hamiltonian gap, 
Aab = (ea|A|eb), and jca), a = {0,1} denotes the instan¬ 
taneons eigenvectors of H{s). 

We now illustrate the QAT for H{s) = uJx{^ — s)a^ + 
LUzSa^ and system operator A = where 5 is a 

coupling constant. For A = ga"^, there is always a 
gap in the Lionvillian spectrum above the zero eigen¬ 
value. Correspondingly the decay is T~^ as shown in 
Fig. 2(a). For A = gcr^ the system Hamiltonian com- 
mntes with A when s = s* = 1, at which point the spec¬ 
trum becomes degenerate. Indeed one can check that 
jAoij^ = a;^(l — s)^/||m(s)p and, correspondingly, from 
Eq. ( 8 ), the second eigenvalue goes to zero quadratically 
ins —s*: A 2 (s) ~ — 7 ( 5 )(a;a;/w^)^(l — s)^(l-be”^"^). This 
corresponds to a = 2 and hence [recall Eq. (5)[ an expo¬ 
nent 77 = 1/3. Onr nnmerical simulations are in agree¬ 
ment with this prediction, as seen in Fig. 2(b). Note that 
the Hamiltonian gap 5 enters only indirectly via A. 


VIII. CONCLUSIONS 

Using a novel adiabatic expansion we have extended 
the adiabatic theorem of quantnm mechanics to open 
systems described by a time-dependent master eqnation 
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with generator C{t) in Lindblad form. The theorem, first 
proven using different methods in [20-22], states that if 
one initializes the system in Ker£(0) and a gap condi¬ 
tion is satisfied, the evolution brings the system close 
to Ker£(T) up to an error C/T, where C is a constant 
and T the total time. Our approach allowed us to ex¬ 
tend the results of [20-22] in two directions particularly 
relevant for quantum state preparation and quantum an¬ 
nealing in open systems. On the one hand, we related 
the constant C to the smallest (in absolute value) gap 
Amin of the Liouvillian. For general Liouvillians we ob¬ 
tained C = 0(A“?j), whereas for thermal baths satisfy¬ 
ing the KMS condition we found an improved scaling 
C = 0(Am?j^). More precisely, we showed that tak¬ 
ing T > /3|i£'(cr)||maxA/([-ff'(f^)]^)G.maxAmL/e, guaran¬ 
tees adiabaticity up to an error 0(e) in trace norm. On 
the other hand, we extended previous results to the case 
of level crossing, for which the error becomes 0{T~^) 
with an exponent ry that depends on the rate of the 
gap closing. Thus level crossings with the instantaneous 
steady state can slow convergence down. We provided 
several examples to illustrate our findings, which con¬ 
firm the predicted scaling with T. An interesting open 
question is whether the growing body of techniques de¬ 
veloped for bath engineering [41-43] can be used to enact 
boundary cancellation methods and reduce the error to 
O (£“”) with controllable n > 1, as in the closed sys¬ 
tem case [11, 44-46]. Our results have implications for 
adiabatic quantum computation and quantum annealing 
in the presence of dissipation, where the closed-system 
adiabatic theorem cannot be directly applied. 

ACKNOWLEDGMENTS 

This work was supported under under ARO MURI 
Grant No. W911NF-11-1-0268 and ARO grant number 
W911NF-12-1-0523. 


Appendix A: Proof that zero is a semi-simple 
eigenvalue of £(s) 

We assume in this section that £(s) is finite¬ 
dimensional. We show that zero is a semi-simple eigen¬ 
value of £(s), i.e., that there are no idempotents in the 
Jordan block decomposition of £ corresponding to its 
zero eigenvalue. Assume the contrary, i.e., that one can 
write £(s) = 0 x P{s) + D{s) + R{s) where D{s) is nilpo- 
tent, i.e., D{s)'^ = 0, £)(s)'"“^ ^ 0, and R{s)P{s) = 
P(s)i?(s) = 0. Then 

= P{s) + (Al) 

n—1 

where we also used that D{s)P{s) = P{s)D{s) = D{s). 
By assumption, £(s) is Markovian, which means that for 
each s and t > 0, is CPTP and so |j = 1. 


Taking the norm of both sides of Eq. (Al) we obtain 
that the left hand side is bounded while the right hand 
side grows unboundedly with t, which is a contradiction. 
Hence D{s) — 0 and £(s)P(s) = P(s)£(s) = 0. For an 
alternative proof see, e.g.. Ref. [47] or Proposition 5 of 
Ref. [22] . 

Note that when the generator £(s) is normal, (i.e., 
[£(s),£*(s)] = 0, where * denotes the adjoint with re¬ 
spect to the scalar product {X,Y)a ■= Tr[yOG(s)-^^y])) 
as is the case for a thermal bath, this has the pleas¬ 
ant consequence that all the eigenprojectors of £(s) are 
bounded (in fact have norm one) and there are no nilpo- 
tent terms in its Jordan block decomposition. 


Appendix B: Proof that P{s)V{s) = V{s)P{0) 

Recall that we defined the intertwiner to be the so¬ 
lution of the following differential equation: P'(s) = 
[P'(s),P(s)]R(s). To check that R(s) satisfies the in¬ 
tertwining relation, define W{s) = P{s)V{s). From 
now on we drop the s dependence when not explicitly 
needed. Differentiating P^ = P gives P' = P'P + PP', 
and after right-multiplying by P gives PP' P — 0 and 
likewise QQ'Q = 0 [which also implies QP'Q — 0 
since 0 = QVQ = Q{P' + Q')Q\. Thus [[P',P],P] = 
P'P — 2PP'P + PP' = P'. Using this, note that W = 
P'V + PV = {P' + P[P',P])V = [P',P]PV = [P',P]W, 
which is the same differential equation as the one satis¬ 
fied by V. Since W and V differ in their initial condi¬ 
tion, i.e., U(0) = 1 and W{0) = P(0), it follows that 
W(s) = V(s)P(0), i.e., 

W(s) = P(s)V(s) = V(s)P(0) . (Bl) 


Appendix C: Complete positivity of U(s)P(0), and 
lack of positivity of U(s) 

Here we show that W(s) = V(s)P(0) is a CPTP map 
even though the intertwiner V(s) itself is not in general. 
From Appendix B we obtain 

W' = [P',P]W 

= P'PW - PP'W = P'PW - PP'PW 
= P'PW = P'W . (Cl) 

Now we write the solution of the differential equation 
using the method of Euler lines [48, 49], i.e., W{s) = 
lim 7 v_,,oo (1 + eP'(s — e)) • • • (1 -I- eP'(O)) P(0), with e = 
s/N. Using P(s)^ = P{s) we get, in the e —>■ 0 limit, 

(1 -f eP'(O)) P(0) = (P(0) -k eP'(O)) P(0) = P(e)P(0) , 

where the second equality is up to first order in e. Then 
(with the same notation) 

(1 + eP'(e)) P(e) = (P(e) + eP'(e)) P(e) = P(2e)P(e) , 
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etc., until W{s) = limAr_>oo-P(s) ■ ■ • P(2e)P(e)P(0). In 
this form, W (s) is an infinite product of CPTP maps, so 
it is a CPTP map itself. It follows from submultiplica- 
tivity of the norm [50] that ||IT(s)|| = ||P(s)P(0)|| < 1. 

Next, let us demonstrate that P(s), in general, is not 
even positive. Consider, e.g., the case where P(s) is one¬ 
dimensional, so that P(s) = |/9(s))(l|, where {x\ is the 
adjoint of \x) with respect to the Hilbert-Schmidt scalar 
product {A\B) = Tr[Alp]. We assume that p{s) is differ¬ 
entiable. Then P'(s) = |/9'(s))(l| and so P'P = P' while 
PP' = 0 since p'{s) is traceless. Thus V' = [P',P]V 
becomes V = P'V (with P(0) = 1), with the solution 

V{s) = Texp P'{a)da^ , (C2) 

where Texp denotes the time-ordered exponential. How¬ 
ever P'(si)P'(s 2 ) = P'{s 2 )P'{si) = 0 [again, since p'{s) 
is traceless] so the time-ordered exponential in the above 
equation reduces to a standard exponential. The solution 
is then 


P(s) = exp[(|p(s)) - |p( 0 )))(l|] 

= l + (|p(s))-|p(0)))(l| , (C3) 

where in the second line we used the fact that [(|p(s)) — 
|p(0)))(l|]^ = 0. Now, it is easy to see that, unless p{s) = 
p{0) such that P(s) = 1 , P(s) does not even preserve 
positivity. Let us define, for clarity, 6p := p{s) — p(0). 
This is a hermitian and traceless operator, so it must 
have a negative eigenvalue, i.e., there exist a vector |a) 
and a real number a > 0, such that 6p\a) = —a\a). Let 
us consider the following state 

a;o := -h (1 - A)|a-^)(a-^|, (C4) 

where d is the dimension of the ISS, is a vector 

satisfying (aja-*-) = 0 and A € [0,1]. It is clear that xq is 
a positive operator. However, 


P(s)a:o|a) = Q - |a) , (C5) 

such that, for A < min(a(i, 1), V[s)xq has a negative 
eigenvalue, i.e., P(s) does not preserve positivity. 


Appendix D: Proof of Eq. (2) 

We assume that C is m-fold differentiable. 


1. The m = 1 case 

The proof of Eq. (2) of the main text for m = 1 closely 
follows the classic reference [4[. As above, let P(s) be the 
projector onto the instantaneous, zero eigenvalue of C{s). 
In this subsection we do not require the ISS to be unique. 


Since we established in Appendix A that CP = PC = 0 
(from now on we omit the explicit dependence on s if not 
strictly needed), it follows that C — QC = CQ. Therefore 
the reduced resolvent S = lim^-^o Q{C — z)~^Q satisfies 
SC = \iva.z^oQ{C-z)~'^QC = Q\\vclz^q{C- z)~^C = Q. 
To summarize: 

CS = SC = Q , (Dla) 

S-P = PS- = 0 . (Dlb) 

It follows from Eq. (Cl) and PP'P = 0 that PW' = 0 
and so 

W' = QW' . (D2) 

Recall [Eq. (1) of the main text] that the evolution 
operator £ satisfies £' = TC£. Now, 0 = {££~^y = 
£'£~^ P £[£~^)y so that 

(£■-1)' = -T£-^C . (D3) 

Therefore {£-^yW = -T£-^CPV = 0, since CP = 0. 
Hence {£~^Wy = £~^W', which integrated over [0, s] 
gives 

£-\s)V{s)P{0)-P{0)= [ £-\a)W'{a)da . (D4) 

Jo 

Now we use 

£-^W' = £-^QW' 

= £-^CSW' 

=-T-\£-^ySW' , (D5) 

where we used Eqs. (D2), (Dla), and (D3). We plug this 
into Eq. (D4) and integrate by parts, to obtain 

£-'^{s)W{s) - P{0) = 

^ £-^{a)[S{a)W'ia)]'da - (^ 6 ) 

Now we act with £{s) from the left and use the property 
£ (a, h)£ {b, 0) = £ (a, 0 ) to obtain 

[£-(s) - E(s)]P(0) = ^ £{s,a)S{a)Wya)\; 

[ da£{s,a)iS'W'+ SW")ia) . (D7) 

T Jo 

At this point recall that W' = P'W [Eq. (Cl)[ and 
hence [using Eq. (Bl)[ W'{s) = P'(s)V(s)P(0), and also 
W" = [P" + {P'Y]W. It follows from the latter, to¬ 
gether with S = SQ [from Eq. (Dlb)[ and W = PW 
[from Eq. (Bl)[, that SW" = SQ[P" + {P'f]PW. But 
Q[P'YP = QP'{P p Q)P'P = 0 because PP’P = 0 and 
QP'Q = 0 (see Appendix B), so that SW” = SP”W. 
Collecting all these results we finally obtain 


[£:(s) - E(s)]P(0) = 

iS{s)P’{s)Vis)Pi0) - £{s)Si0)P’{0)Pi0)) 



da£{s,a) [S'P' P SP”] (tT)E(a)P(0)} . 


(D 8 ) 
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Using submultiplicativity along with ||P(s)|j = 1 and 
|jU(s)P(0)|| < 1 (as shown in Appendix C) and ||£(s)|| < 
1 (see immediately below) we obtain 

\\[£{s)-V{s)]Pm<C/T , (D9) 

with 

c = |!5(s)|!||p'(s)|| + ||5(o)||||p'(o)|| 

+ sup \\[S'P' + SP"]{a)\\ , (DIO) 

<tG [ 0 , s ] 


as stated in Eqs. (3) and (4)of the main text. 

Note that when restricted to acting on normalized 
states, clearly Ill’ll = 1 by trace preservation. More gen¬ 
erally, since C{s) is a Lindbladian we have |j exp[t£(s)] || < 
1 for each fixed s and t > 0 [i.e., C{s) generates a con¬ 
traction for any fixed s]. Using the method of Euler lines 
(see Appendix C) one can show that this implies that 
ll'^^(sU)ll ^ 1 ) because £{s,t) can be written as an infi¬ 
nite product of evolutions of the form exp[t>C(si)]. 


2. Extension to arbitrary order 

Next we show how to extend the integration by parts 
technique of the previous subsection to arbitrary order 
m. A particularly simple series will result under the as¬ 
sumption that the ISS is unique, i.e., P{s) = |p(s))(l|, 
as in Appendix C (this assumption is not essential, but it 
significantly simplifies our calculations below). As before, 
we also assume that C generates a trace preserving map. 
Both of the latter two assumptions are usually satisfied 
for thermal (i.e., Davies) generators. 

We begin by considering the first term of Eq. (D6), i.e., 
{SW'Y = S"W' -f SW”. Our strategy will be to repeat 
the integration by parts on the term arising from W', 
while keeping (i.e., not integrating by parts) the term 
arising from W”. We will show how this can be done 
repeatedly and thus extend the result to arbitrary order. 

Now, thanks to the assumption of uniqueness P{s) = 
|p(s))(l|. Hence from PS = 0 we obtain (IIS' = 0 and 
so P'S = |p'(s))(l|S = 0. Using P'S + PS' = Q we 
get PS' = 0 so that, writing S'Q = {P + Q)S'Q = 
PS'Q + QS'Q it thus follows that 

S'Q = QS'Q . (Dll) 

This means that, using Eq. (D2), we can write 


where in the last equality we used Eq. (D3). This term 
can now be integrated by parts again, in analogy to the 
term in Eq. (D5). Differentiating SS'W', we obtain two 
terms, one of which is SS'W", which we keep. The other 
is {SS'yW' which we would like to integrate by parts 
again (because it diverges more strongly when Amin —>■ 
0). We now show that this procedure can be iterated to 
any order provided C is differentiable sufficiently many 
times. 

Let us assume that at order n we obtained the term 

^ l^^da£-^X'„W'. (D14) 

Assume for the moment that = QX'^W'; we 

will show shortly that this is legitimate. In this case 
we can repeat the calculation of Eq. (D13) and write 
£-^X'^W' = -T-^(£-Q'SX'^W'. We can thus use the 
integration by parts trick: 


— / du£-^X'J¥' = - 


/o 


1 


1 


?-l CV' w'l 


£-^sx'„w'\: + 


Tu+i 


1 


2^n+l 

+ 

1 


f da£-^{SX'„W')' 
Jo 


da£-^SX'„W" 


^ lo ' 2"n+l 




(D15) 


The first integral (with lU") is the one we keep, while the 
second (with W') is the one we continue to integrate by 
parts. This shows that the process can be iterated with 
A„+i = SX'„. 

We now verify that we can plug in a Q term at each 
order as claimed above. We prove it by induction. As¬ 
sume that X'^W' = QX'^W'. We wish to show that this 
implies that X'^j^^W' = QX'^j^iW'. But this is clear 
since 


= {s'x'^ + sx'Y)w' 

= {S'QX'^ + QSX'Y)QW' = Q{S'QX'^ + SX'Y)QW' 

= Q{S'X'^ + SX'Y)W' = QX'^^.W' , (D16) 

where we used the induction hypothesis, Eq. (D2) and 
Eq. (Dll). We have already shown that the claim holds 
for n = 1 [Eq. (D12) for Xi = so we are done. 

We thus obtain 


S'W' = S'QW' = QS'W' . (D12) 

Returning to the term of interest in Eq. (D 6 ), we have 

£-^S'W' = £-^QS'W' (D13a) 

= £-^CSS'W' (D13b) 

= -T-\£-^ySS'W' , (D13c) 


^ -I nS 

e-'(,}w{s) - pm = i: A + ™ 

n=l ^ ^ Jo 

(D17) 

Tr, = -£-^X^W'\lp [ da£-^X^W" 

(D18) 

A„+i = SX'^ , X^ = S. (D19) 








Multiplying from the left by £{s) we obtain 


Appendix F: Extension to the case of a degenerate 
kernel 


[£(s) - V{s)]P{0) 




m ^ 

E i Lfi 

'T'n 


^ j da£{s,a)X'^{a)W''{a) 

(D20) 

£:(s,cr)A„(cr)VF'(cr)|o 

[ da£{s,<T)Xr,{a)W"ia). (D21) 
Jo 


If the kernel of C{s) is not one-dimensional, in general 
we will also have a P-Q block in The procedure, 
however, can still be iterated. Assume that at order n 
we obtained 

da£-^X„W' = 

da [£-'^PXnQW' + £-^QXr,QW'] . (FI) 



If C is differentiable infinitely many times, we can write 

°° O 

[£-(s) - V{s)]P{Q) = ^ ^ , (D22) 

71 — 1 

under the assumption of convergence. We expect this 
assumption to hold in the finite-dimensional case if C 
depends analytically on s. 


Appendix E: Additional remarks on the closed 
system limit 


The closed system limit may also be achieved by con¬ 
sidering £ 3 ,( 5 ) = K,{s)-\-x'D{s), where I?(s) is a dissipative 
generator, and taking the (weak coupling) limit a: —>■ 0 . 
This is a singular limit: for x ^ Q there is typically only 
one steady state, whereas for a; = 0 there is a large degen¬ 
eracy (> d Hilbert’s space dimension). Calling P{x) the 
spectral projection onto the zero eigenvalue for fixed s, 
one has lima;_,.o P{x) yf P{0), since the left hand side has 
rank one, whereas P(0) has rank > d. Moreover, for suf¬ 
ficiently small X the smallest gap will be oc a; and hence 
also (1/gap)’* {rj > 0) will diverge. Thus the constant C 
in Eq. (.), and even the optimal constant Copt [obtained 
taking the infimum of all the constants C satisfying the 
bound Eq. (o) for all T] will, in general, likely be large. 
In other words, in general the limit a; —> 0 is singular 
and the adiabatic theorem for closed system cannot be 
recovered in this form. A possible way to recover it is to 
abandon projectors and move to states, at least in the 
case of non-degenerate ISS for a; yf 0. In other words, 
let us call Px{s) the state evolved with generator Cx{s) 
from ax{0) [£x(s)cr 2 ,(s) = 0 is assumed to be unique for 
a: yf 0]. Now it is possible that Px{s) and ax{s) have a 
well defined limit as a: —)■ 0 such that the following bound 

|jp.(s)-a,(s)|| <a/T (El) 

may admit a non-trivial limit for a: —>■ 0 with 

lima._,,o Cx < 00 . This may be achieved by applying 
Eq. (D 8 ) to ax{0) and bounding the resulting expression. 


We keep the first term (we will have to show later 

that it scales nicely), and integrate the second by parts, 

for which we use £~^QX„QW' = £~^£SX„QW' = 
—T~^(£~^)'SXnW'. The above equation becomes 

+ ^ iya£-\SX„W'y (F2) 

+ ^ ly a £-\SX„)W'' 

+ iya£-\SXjW' , (F3) 

from which we see that the recurrence is A„+i = (SX^)' 
with Xq = 1 as previously. However, now the series reads 


00 -p, 

£-\s)W{s)-P{0) = J2^ (F4) 

n—1 

rn = -£-^x„w'\l 

+ [ da£-^XnW" 

Jo 

+ [ da£-^PXyQW' , (F5) 

with Xi = S. Multiplying from the 

left by we obtain 




[£(s) - I/(s)] P(0) 

n—1 

£ln = £{s, a)Xnia)W'{a)\l + 

— I da £{s, a)Xn(a)W"(a) 

Jo 

- r da £{s,a)[PXyQW']{a) . 
Jo 


(F 6 ) 


(F7) 


In order to repeat the arguments for the case of level 
crossings, we now have to assess the scaling of PX'^Q as 
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one eigenvalue goes to zero. We have separately shown 
that Xn 1/(5^’* with /?„ = na + n — 1. Now, = 
5(x;_i), so 

PX;Q = P(5(x;_i))'Q (F8) 

= PS'X'„_,Q (F9) 

= PS'QX'^_,Q , (FIO) 


where the first equation holds because PS = 0 and the 
second because PS'P = 0. Now, as we know, PS'Q = 
—P'S, introduces only a power 1/(5“. On the other hand 
we know that X'^_i ^ Hence 


PKQ‘ 


1 

(5/^n. —1 + 1+0 




(Fll) 


We see hence that this potentially more dangerous term 
has the same scaling as Xn- In other words, all our ear¬ 
lier conclusions also hold for the case with a degenerate 
kernel. 


Appendix G: Proof of an identity for S' 

Since S = SQ [Eq. (Dlb)] we have S' = S'Q + SQ'. 
Note that 

SQ' = -SP' = -SQP' = -S^CP' = S^C'P , (Gl) 

where we used that SC = Q, and CP' = —C'P. Differ¬ 
entiating SC — Q yields S'C + SC = —P'. Multiplying 
from the right by S we obtain 

S'Q = -P'S - SC'S . (G2) 

Therefore we have: 

S'Q = -P'QS - SC'S 

= -P'CS^ - SC'S = PC'S^ - SC'S . (G3) 

Gombining the last three equations, we have: 

S' = S^C'P + PC'S^ - SC'S . (G4) 


Appendix H: Gap dependence for unitary families 

Assume here that C{s) is a unitary family, i.e., C{s) = 
with /C an anti-hermitian superoperator. In 
this case one has (dropping the s dependence) C = XC — 
CX and so, using first order perturbation theory ([35], 
page 77): 


P' = -PC'S - SC'P (HI) 

= -PXQ + QXP (H2) 

= XP-PX . (H3) 

Differentiating the above equation we get that both P' 
and P" are sums of products of P and X and so are 
everywhere bounded under the assumption that P and 


X are bounded. From Eq. (G4) and repeatedly using 
SC = CS = Q, CP = PC = 0, and S = SQ = SQ [from 
Eq. (Dlb)[ we obtain 

S' = S^{XC - CX)P -f P{XC - CX)S^ - S{XC - CX)S 
= -SQXP + PXQS - SXQ + QXS 
= -SXP + PXS - SXQ + QXS 
= -SX + XS . (H4) 

Hence in this case S' carries only one inverse Liouvillian 
gap (hidden in S). Gombining these results we see that 
the constant C appearing in Eq. (DIO) is in this case at 
most 0(A”(„). 


Appendix I: Bound on | P/ for the KMS case 


In the KMS case with a unique steady state we have 
already seen [Appendix C[ that P{s)x = pc{s)Tt{x) with 
Pq = /Z and the partition function Z = . 

Then obviously P'{s)x = p'q{s)Tt{x) and ||P'|| = ||/9'||i. 
We first recall the Bogoliubov-Duhammel scalar product 
[51, 52[: 


(++ = 


dx Tr 




(II) 


and define {A)g ■= Tr [Ae jZ. Moreover one has 
the inequality [52[ 

(A,A) < i(AlA-f AAl) . (12) 

Differentiating the Gibbs state we obtain 

Pg = dx + ^^pg{H') , (13) 

where the second term arises from the differentiation of 
the partition function. Now let Y denote the first term 
of the above equation. We wish to bound its trace norm. 
We use [50[ 

liFIli =sup|Tr(yC/)| , (14) 

u 

where the supremum is taken over the set of unitary ma¬ 
trices U. Note that Ti (YU) = —P{H',U), since H' 
is Hermitian. Using the Schwartz inequality we have 
\Tr{YU)\ < P^{H',H') {U,U). We now use Eq. (12) 
to obtain 


nil <. 


(15) 


On the other ||/0g|| = 1- Gombining the bounds we get, 
all in all 


iip'iii<^(+((nn + i(nG 


(16) 


Since in general > (A)^, we can write an even more 

compact bound : 


iip'ii = iip'iii<2/3v+nn 

as reported in the main text. 


( 17 ) 
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Appendix J: Extension to the case of level crossing 

The results obtained so far assume that there is a fi¬ 
nite gap from the zero eigenvalue to the rest of the spec¬ 
trum. However, this is not necessary and all of the above 
can be generalized to include the case in which there 
are level crossings along the adiabatic path [0, s]. Sup¬ 
pose that there are N level crossings at positions s*, 
i = 1,..., TV. We isolate each level crossing by a short 
segment of length 5i, so that no two segments overlap. 
Then define two sets A, B such that [0, s] = A U H, with 
^ = Ui['Si+<^*/2,s*+i-<5i+i/2], i = 0,... ,TV, with = 0, 
s^+i = s, (5o = 5n+i = 0, and B = \Ji[s*-5i/2, s*+5il2], 
i = l,...,iV. By construction, A avoids the singulari¬ 
ties while B comprises short segments around each level 
crossing. We start with Eq. (D4), which we separate into 
integrals over A and B: 


s*)“' (where the coefficients Vi can be complex). This 
algebraic vanishing of the gap is the only possibility in 
the finite-dimensional case on which we focus. Moreover, 
for normal operators and analytic dependence on s, the 
exponents ai must be positive integers (see Ref. [35]), 
but we won’t be needing this. Consider then the finite¬ 
dimensional case. The normality of C guarantees that 
there are no nilpotent terms in the spectral resolution of 
C and so the reduced resolvent can be written as 

s(») = E yi ■ 

i>o ^^ ’ 


where Aj(s) [Rj^s)] are the eigenvalues [spectral projec¬ 
tions] of £(s), with Po(s) = P(s). Hence 


j>0 ^ ' j >0 


(J6) 


£-\s)V{s)P{0)-P{0)= [ £-'^W'da 

Jaub 

= -^ [ i£-^ysw'da+ [ £-^W'da , (Jl) 
^ J A J B 

where in the integral over A we used Eq. (D5). From the 
above equation, using the classic argument in Ref. [4] one 
can than show that 

lim [£(s)-R(s)]P(0) =0. (J2) 

T—>00 

However, we are interested in assessing the way in which 
the limit is approached. In order to obtain a better es¬ 
timate of the asymptotic behavior it is necessary to re¬ 
peat the integration by parts trick at all orders. Since in 
region A there are no exceptional points we can imme¬ 
diately repeat all the steps from Appendix D. Assuming 
convergence of the series we then obtain 


Now, it is natural to assume that Pj are piecewise dif¬ 
ferentiable also for j > 0. For normal operators and 
analytic dependence on s this follows from a Theorem of 
Kato (see Ref. [35] Theorem 1.10, page 71). In this case 
the largest divergence of S' , when s — s* = —>■ 0, comes 
from the vanishing of the gap in the first term of Eq. (.16) 
and is of the form 

S'(») ~ jy+i ■ (J7) 

Note that the prefactors may be different depending on 
whether 5i -A 0^. This reasoning can be extended to 
A„(s). Recall that A„ = with Xi = S, and 

note that A„ contains n powers of S and a total of 
n — 1 derivatives with respect to s. So we conclude that 
-^n(s) i^^nui+n-i gg , 5 ^ More precisely we have 

Jim df-XUs* + S,)=Y^ , (J8) 

Oi —>-0^ 


qA n 

[£[s)-V{s)]P{0) = Y,^- / da£{s,a)W'{a) 

n=l ^ ^ 


N 


i=0 


st+i-Si+i/2 

s‘+5i/2 


H)^ = ^£(s,a)A„(n)lE'(a) 

) 

— J da £{s, a)Xn{cr)W" (a) . 


(J3) 


(J4) 


where the Yj exist and where 

= nai + n-l . (J9) 

Thanks to the above there exist positive constants 
such that, for sufficiently small di 

A, A* 

i=\ °i " 


Each term in the above equations is bounded for finite 
5i, however, when (5^ —> 0 some terms display singular 
behavior. We are interested in estimating the size of 
the most singular terms. Recall that we are assuming 
a Davies generator. In this case W = P'W [Eq. (Cl)[ 
is bounded (because P' is bounded [32] and ||1T|| = 1 
because W is a CPTP map), and similarly W" = [P" + 
{P'J]W is bounded because P" is bounded [32]. Hence, 
when 5i -A 0, the only singularities in Eqs. (.I3)-(^'') 
arise from A„((t). Now assume that at each singular 
point s*, there is an eigenvalue going to zero as Vi{s — 


The A\^ can be made independent of 5i, for sufficiently 
small 5i. In fact they can be taken to be the max of the 
norm of the most diverging terms plus a small constant. 
On the other hand, to bound the second term of Eq. (.13) 
we just need to notice that W is piecewise continuous 

II - f da£{s,a)W'{a)\\ 

J B 

= ^ da£{s,a)W'{a)\\<Y,^^B,, (Jll) 

i=l i=l 
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with positive Bi, because the integrand is piecewise con¬ 
tinuous {B, - (l/2)(||P'«+0+)|| + ||P'(s*-0+)||). All 
in all we obtained 


N 

|l[£(s)_y(s)]P(0)|! <^/,(r,50 , (J12) 

2=1 


where 


n—1 



+ Bi5i . 


(J13) 


obtain 


hiT,6^in) = 


D, = 


D, 


{c^TY 




.n—1 


Returning to Eq. (.112) we finally obtain 


N 


[£{s)-v{s)]pm<j2 


D, 


as stated in Eq. (5) of the main text. 


(J16) 

(J17) 


(J18) 


Now for each fixed T, we are free to minimize each fi 
over 5i, to which we turn next. 

Differentiating, we obtain 


dsji 


E 




-\- Bi — Q , 


(J14) 


Appendix K: Derivation of the eigenvalues of the 
Lindbladian of Example 2 

Denote the instantaneous Hamiltonian gap by (5(s) and 
the instantaneous ground state and first excited state by 
|eo(s)) and |ei(s)) respectively. The action of the Lind¬ 
bladian on p can be written in this basis as (dropping the 
s dependence for clarity): 


where /3(i + 1 = n{ai + l). Recall that all the > 0 and 
P* > 0, so the only way Eq. (.114) can have a solution is 
that, asymptotically as T —>■ oo, some of the terms in the 
series in Eq. (Jl-i) scale to constant, while some others 
may scale to zero. Plugging in the Ansatz Si = {ciT)~'^' 
one can check that pi = 1/(0^ -I- 1) is the only possible 
solution. In fact all the terms scale to a constant. The 
constants Ci must satisfy 


Y,PnAlc^ = B,. (J15) 

n=l 

The extremum is indeed a minimum as it is easy to check 
that the second derivative is always positive for positive 
T,S. Substituting the solution back into Eq. (.113) we 


{eo\Cp\eo) = "/{S) [|Aoippii - e"^^|Aioppoo] (Kla) 
{eo\Cp\ei) = (ip-T) poi (Klb) 

(ei \Cp\eo) = {-ip - E) pio (Klc) 

{ei\Cp\€i) = j(S) [e"^^|Aioppoo - |Aoippii] , (Kid) 

where pab = {ta\pYb), p = 6 - S'(0)(A§o - + 

[5(5) - 5(-5)]|Aoi|2 and 2r = 7(0)(|Aoo|" + |An|") + 
7 ( 5 ) |Aoi|^ (1 + This form immediately allows us 

to read off two of the eigenvalues of £ as zLip — T. The 
remaining two are the eigenvalues of the following 2x2 
matrix: 

7(5)|AioP , (K2) 

whose eigenvalues are 0 and —|Aoip7(5) (l -I- 
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